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Abstract: We use perturbation theory to compute the diffusion coefficient of a heavy 
q \ quark or scalar moving in M = 4 SU(A r c ) Super Yang-Mills plasma to leading order in 

the coupling and the ratio T/M <^ 1. The result is compared both to recent strong 
coupling calculations in the same theory and to the corresponding weak coupling result 
in QCD. Finally, we present a compact and simple formulation of the Lagrangian of our 
o theory, M = 4 SYM coupled to a massive fundamental Af = 2 hyp er mult iplet, which 



is well-suited for weak coupling expansions. 
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1. Introduction 

An important and interesting challenge facing theorists investigating heavy ion physics 
is to predict the rate of energy loss of a heavy quark moving through a quark gluon 
plasma (QGP), which is a quantity of direct experimental relevance Jl|, |2], [3]. For weak 
coupling and ultrarelativistic quarks, 'jv ^> 1/g ^> 1, the dominant mechanism for this 
is gluon bremsstrahlung, while for the experimentally equally important region ^yv < 1, 
the energy loss occurs through elastic collisions with light plasma constituents. Both 
of these cases have been extensively studied || 0, || , but all existing calculations 
share the fundamental shortcoming that they assume the plasma to be weakly coupled, 
which need not be the case in the energy range of interest for instance for RHIC physics. 

While one has grown to rely upon lattice QCD as the source of direct information on 
the strong coupling regime of various static observables, it is at present still a relatively 
inefficient tool in the description of real-time phenomena (for recent advances, see e.g. 
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Ref. ||). A lot of attention has therefore been turned towards addressing the question 
of the energy loss rate of a quark moving in a strongly coupled non-Abelian plasma in 
an entirely different framework, the M = 4 supersymmetric Yang-Mills (SYM) theory 
with gauge group SU(iV c ). There, one has the unique opportunity of being able to 
access analytically the strong coupling limit of the quantum field theory (in its large 
N c limit) via the famous AdS/CFT conjecture that relates the theory to dual type IIB 
supergravity on an AdS$ x S 5 background [JTOj . The energy loss calculation then reduces 
to studying classical string dynamics in the AdS^ background, which has yielded useful 
analytical results for the heavy (and light) quark energy loss parameters in the strong 
coupling limit (see Refs. fTl], [12], [13], |14| [15], |ll| |l7j and references therein). 



While the development of non-equilibrium lattice field theory methods as well as 
the search for dual string theories of more QCD-like theories continue, it is worthwhile 
to first ask the more modest question of what kind of qualitative, or even quantitative, 
insight the QCD community can draw from the existing SYM calculations. An obvious 
way of addressing this is to perform similar weak coupling calculations in the SYM 
theory that have been carried out in the QCD context and compare the results on 
one hand to the the weak coupling limit of QCD and on the other hand to the strong 
coupling limit of the M = 4 SYM. This can be expected to yield valuable information 
on the similarities and differences of the two theories and to furthermore indicate, to 
which extent one can simply extrapolate the existing weak coupling results in QCD to 
strong coupling. In the present paper, we aim to do exactly this by investigating the 
simplest observable related to the energy loss of a non-relativistic heavy quark in the 
SYM theory, its diffusion coefficient, in the weak coupling regime. The calculation is 
to a large extent parallel to the corresponding QCD computation of Moore and Teaney 
H and generalizes many of its results to the SYM case. 

Our paper is organized as follows. In Section 2, we briefly introduce M = 4 SYM 
and write down its Lagrangian in a form useful for weak coupling expansions, after 
which we discuss when and how one may use semi-classical kinetic theory techniques 
to obtain the diffusion coefficient of a heavy particle in a quantum field theory. In 
Section 3, we then go through the necessary calculations, after which we display our 
main result for the heavy quark diffusion coefficient, Dq, that is to be further analyzed 
and discussed in Section 4. In Section 5, we finally draw conclusions and outline 
some future work to be carried out through weak coupling kinetic theory calculations 
in the SYM theory. The Appendices contain further computational details, such as 
a relatively detailed derivation of our Lagrangian as well as a discussion of how the 
necessary scattering amplitudes and integrals were evaluated. 

Throughout the paper, we will for convenience use notation adopted from Refs. M, 



IS]. This implies working with the Minkowski space metric h ++ and denoting 
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four- vectors by capital letters P, Q, three-vectors by bold ones p, q and the absolute 
values of the latter by p,q. The Dirac gamma matrices are defined so that 70 is anti- 
hermitian while the 7, are hermitian, and consistently with this we have 75 = ^7o7i7273 
and (3 = ijq. The gauge is fixed by choosing to work in the Coulomb gauge. 

For the most part, we will work with four-component Majorana spinors that can 
be written in the special form 



c 

where e = 102 and ( denotes a two-component Weyl fermion. For the Majorana spinors 

$ = tfP = ^ T e 75 , (1.2; 

with e being the 4x4 matrix 

fe 
6= U e 

2. The setup 

2.1 J\f = 4 SYM with massive quarks 

The theory we consider is M = 4 Super Yang-Mills coupled to a M = 2 heavy funda- 
mental hypermultiplet. The field content of M = 4 SYM consists of a gauge field A^, 
four Majorana fermions ipi and three complex scalars <ft p , while the additional M = 2 
multiplet is composed of two heavy scalars $ n and a Dirac fermion u. All M = 4 
fields transform under the adjoint representation of the gauge group SU(iV c ) and are 
therefore traceless hermitian N c x N c matrices, while the M = 2 sector consists of N c 
component vectors in color space transforming under the fundamental representation. 
Following Ref. |L9|, we define <p p = 1 / \/2 (X p + iY p ) , with X p and Y p hermitian, 



which allows us to write our Lagrangian in the form 

C = C + C 1 + C 2 , (2.1) 



with 



C = -tr + ^Mi + {DX p f + (DY P ) 2 } 

- & n (-D 2 + M 2 )$ n -O(J0 + M)u, (2.2) 
£1/9 = tr { - i$i a ij [ X p, ^j) + Alsffij^p, i>j] \-&(Yi- 275^1) w 

+ 2V2 Im( - (I>P + ^i$i - $5V>iP + a; + $\4>2P+u} - U)P + i) 2 §2^} 

- 2M& n Y 1 * n , (2.3) 
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C2/9 2 = -\tr (z[ X a,Xb}) 2 + (-l) n $l ([02,4] + lfa,4}) $ n 
- 4Re ($\[<h,<f>3]*2) ~ \ \(-l) n & n t a ®n\ 2 - 2|$ 2 ^i| 2 

(2.4) 

Here, D denotes covariant derivatives in the appropriate representations of SU(iV c ), 
X = (Xi, Yi, X 2 , Y 2 , X 3 , y 3 ) and a sum over repeated indices is implied. The matrices 
a p and f3 p are given by 

^ = (7 ° V « 2 = ( ° 7*1 , « 3 = f ° ?1 , (2-5a) 



icr 2 / \a"i / \— (T 3 

ffl = /-toi \ ;32 = / -ia 2 \ 3 = (- ^ 
V za 2 y \- zcr 2 / V-ctq OJ 

and they satisfy the algebra 

{a p ,a q } = {P p ,p q } = -25 pq , 

[a p ,(3 q ]=0. (2.6) 

For more details on the derivation of this Lagrangian, see Appendix |A|. 

From the form of the Lagrangian it is clear that neither the heavy fermion nor the 
heavy scalar number is independently conserved, as Eq. ( |2.3|) fails to be invariant under 
the separate global U(l) transformations 

uj -> U u u, (2.7) 

where U& and U u are independent phases. This implies that the diffusion coefficients for 
heavy quarks and heavy scalars are in general not independently well-defined. However, 
if U§ = U^, the transformation of Eq. Q2.7p does leave the Lagrangian invariant, which 
means that this combined transformation gives rise to a conserved heavy flavor current 
that includes both fermions and scalars. It is the diffusion of this heavy flavor density 
that we are interested in. 



2.2 Diffusion of a heavy non-relativistic particle 

Following the approach of Ref. f8|, let us consider the kinematics of a heavy particle 
immersed in weakly coupled M = 4 SYM plasma at temperature T. We assume that the 
particle is in thermal equilibrium with the plasma and that its mass M ^ T, in which 
case the typical energy of all types of excitations is E ~ T and the typical momentum of 
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the heavy particle is p ~ V MT. At weak coupling, the dominant scattering processes 
for the heavy particles are 2 <-> 2 elastic collisions with light plasma constituents, in 
which the typical momentum exchanged is q ~ T and the typical change in the heavy 
particle velocity is 5v ~ T/M <C 1. It thus takes many collisions for the velocity 
to change significantly, and consequently the collisions may be treated as uncorrelated 
events, in which the heavy particles receive random kicks from the medium. In addition, 
the mean free path of the heavy particle A MFP ~ (^) g4T | og i is parametrically large 
in comparison with its de Broglie wavelength, thus allowing one to use semiclassical 
methods in the treatment of its dynamics. 

Let us then look at the trajectory of a heavy particle that starts from the origin at 
t — and denote its position at time t by x(£). Under the assumption that collisions 
with the medium force the heavy particle to undergo a random walk, the diffusion 
coefficient D is defined by 

(x 2 (t)) = 6Dt. (2.8) 

Denoting the random force exerted on the particle by the medium by £(t), the above 
assumptions imply that the dynamics of the particle follow from the Langevin equation 
(see e.g. Ref. J2g) 

d ^ = m-m, (2-9) 

where /i is the momentum drag coefficient. As collisions with the light particles are 
uncorrelated, we furthermore have 

(Zi(t)i j ({)) = K6 ij 6(t-1?), (2.10) 

where 3k is the mean squared momentum transfer per unit time. Using the equilibrium 
relation (p 2 ) = 3MT as well as the solution to the above differential equation, it is then 
easy to show that f| 

a = , (x 2 (t)) = —, (2.11) 

/* 2MT V V ;/ M/i' V ; 

from which it follows that the heavy particle diffusion coefficient is given by 

2T 2 

D = —. (2.12) 

K 

The above result relating the diffusion coefficient to k proves highly useful for our 
purposes, as in the semiclassical regime where kinetic theory is valid we may immedi- 
ately write down an expression for the latter in terms of the scattering amplitudes of 
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the quantum theory. Denoting the heavy particles by if, the light particles by £ and 
the Bose and Fermi distribution functions by n^k) and rif(k), respectively, the mean 
squared momentum transfer per unit time is given by @, pTf| 



16(2vr) 5 M 2 



1 




(p - p' + k' - k)5(fc - k') 



x 



H£-^H'v\ 2 n e (k)(l±nt(k')) 



(2.13) 



Hl,H'£' 



Here, I-Mjk-^'^I stands for the scattering amplitudes squared — summed over all 
internal degrees of freedom of the light particles and the final state heavy particle 
and averaged over those of the initial heavy particle (including the flavors) — for the 
process 1 Hi — > H'£' . The plus sign is taken when £' = b represents a final state boson 
and the minus sign when £' = f represents a final state fermion. 

3. Calculations and results 

In the non-relativistic limit where M>T, the number of interaction terms in the La- 
grangian relevant for the scattering of massive particles can be greatly reduced, as there 
exists a hierarchy in the M dependence of the various scattering amplitudes. First of 
all, the amplitude for any tree level process that contains an intermediate heavy parti- 
cle will be suppressed by an inverse power of M relative to those with an intermediate 
light particle. Second, in the non-relativistic limit each external heavy fermion 2 will 
introduce a factor of vM to the amplitudes, and each heavy scalar/gluon vertex (with 
momentum P ~ (M, 0) flowing through it) as well as each heavy scalar/light scalar 
vertex will introduce an explicit factor of M. Therefore, at leading order in g the dia- 
grams proportional to the highest power of M are those, in which a heavy fermion or 
scalar scatters elastically off of a light plasma constituent via the exchange of a light 
intermediate boson. These processes are depicted in Fig. l.a-b, while some examples 
of processes whose amplitudes are suppressed by positive powers of T/M are shown in 
Fig. I.e. The latter include inelastic processes, diagrams with heavy intermediate lines 
and graphs containing a four-scalar vertex that is independent of M. 

Finally, we note that in the non-relativistic limit we may neglect the coupling of 
X p to heavy fermions, which follows from the fact that uj^^uj is parity odd. In the non- 
relativistic limit, the scattering of fermions via scalar exchange is independent of spin 

1 The field theory also allows for heavy particles to scatter off of other heavy particles, but their 
contribution to the integral of k is suppressed by an exponential of M/T . 

2 We use the convention u s (p)?v (p) = 2M5 ss i in the normalization of the heavy spinors. 
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Figure 1: a) The lowest order tree- level processes contributing to the heavy quark diffusion 
coefficient, the amplitudes of which come with the maximal power of M. b) The correspond- 
ing diagrams relevant for the heavy scalar diffusion coefficient, c) Examples of processes 
contributing to the heavy flavor diffusion coefficient but suppressed by powers of T/M. The 
solid and dotted bold lines correspond to the heavy quarks and scalars, respectively, the solid 
and dotted light lines to the corresponding massless fields, and the wavy lines to gluons. The 
arrows drawn adjacent to the light quark lines indicate the (arbitrary) direction of the flow 
of the Majorana fermion number (see e.g. Ref. |22|| ). 



and thus conserves the parity of the fermions, implying that the heavy fermion scatter- 
ing amplitudes containing an exchanged X p must be suppressed by inverse powers of 
M. In summary, to obtain the desired heavy flavor diffusion coefficient to leading order 
in T/M and g, we may neglect £2 entirely and replace Ci by the effective interaction 
Lagrangian 

Cf/g = tr (^75^1, ^]) - - 2M$J,y 1 $ ft . (3.1) 

After neglecting the species changing Yukawa terms from Ci, the Lagrangian be- 
comes invariant under the separate U(l) transformations of Eq. ( |2.7| ). It then follows 
that to leading order in M, the heavy fermion and scalar currents are independently 
conserved, and therefore the corresponding fermion and scalar diffusion coefficients Dq 
and D$ can be independently denned. For these currents, the mean squared momentum 
transfer per unit time is given by 

I f d 3 kd 3 k'd 3 p' 3 



x 

Ll 



{ \Mm-+Ht>\ 2 n t {k){\ ± n e (A/))} , (3.2) 
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where H stands for either Q or S, and the corresponding diffusion coefficients read 



Even though Dq and D$ are independently well defined, a quick inspection of the 
forms of Cq and Cf^ reveals that their values are in fact equal. To see this, note that 
if the momentum exchanged in a collision is q = p' — p ~ T, where p and p' are the 
momenta of the incoming and outgoing heavy quarks, respectively, then the spinors 
u s (p) corresponding to the external legs satisfy (up to 0(T/M) corrections) 



implying that the contributions of the heavy fermions to scattering amplitudes are spin 
independent. To obtain the tree level scattering amplitudes that we are interested in 
(and that are not sensitive to particle statistics), we may therefore replace the heavy 
quark by a complex scalar field E, with the factors of ±2M added explicitly to the 
corresponding vertices. The couplings of E to the light fields are then identical to those 
of $ n , and therefore the corresponding scattering amplitudes agree as well. In what 
follows, we will exploit this symmetry and only consider the heavy fermion diffusion 
coefficient. 

Before we can proceed to the actual computation of Dq, we must still deal with 
the fact that the expression for k q given in Eq. ( |3.2| ) is infrared sensitive, which can 
be seen by noting that if one were to use bare propagators for the exchanged bosons 
in the scattering amplitudes squared, the resulting integrals in Eq. ( |3.2|) would diverge 
in the q — > limit. This can be attributed to the long range potentials associated with 
the exchange of the massless bosons which, however, are modified by the interactions 
with the plasma that cut the divergences off at the scale gT. Taking the interactions 
into account, the IR problem is naturally dealt with by including self energy correc- 
tions to the corresponding propagators. As kinematics furthermore require that the 
energy exchanged in a collision be suppressed relative to the spatial momentum by a 
factor of y/T/M, we note that the appropriate self energy corrections are those due to 
static thermal screening. We can therefore simply add static screening masses to the 
propagators in question, and taking into account the fact that only the temporal gluon 
propagator enters the calculations, obtain as the required resummed propagators 



2T 2 



(3.3) 



M s (p)7 M w s /(p') -2iM8 ss/ 8fj,o, 



(3.4) 



$(P) 



1 




(3.5) 



Gl 



Up) 




p 2 + 171- 

1 



D 



(3.6) 



V 
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for A and <j) p , respectively. 

The squared scattering amplitudes for the processes shown in Fig. |].a are computed 
in Appendix [FJ. Denoting heavy quarks by Q, light fermions, scalars and gluons by /, 
s and g, respectively, and initial and final light particle three-momenta by k and k', 
the results read 

\M Qf 



\M Qs 
\M Qg 
\M Qs 
\M Qg 

where d A = N 2 — 1 and 6 is the angle between k and k'. These expressions have been 
summed over all internal degrees of freedom of the light particles as well as over those 
of the final state heavy quark, and averaged over those of the initial state heavy quark. 

Appendix B shows how to evaluate the integrals appearing in Eq. (|3.2|) . At weak 
coupling, where one may use the leading order results (see e.g. Ref. jl9| ) 

ml = 2g 2 N c T\ (3.8) 
ml = g 2 N c T 2 , (3.9) 

consistency in the weak coupling expansion requires that the terms in the integrals 
proportional to positive powers of m D /T or m s /T be neglected, which allows us to 
carry out the integrations analytically. Using the relation of Eq. (|3.3| ) between k q and 
.Dq, we then obtain as our main result 

12tt f 2T 13 1 C'(2) 1 _1 , 

D Q = -—-\ log— +- - 7E + -log2 + y-i \ , 3.10 
d A g 4 T I m D 12 3 C(2) J 

to which the leading corrections come in at relative order 0(g). This result is indepen- 

2 

dent of the scalar screening mass, which follows from the fact that as cos 6 = 1 — > 
every term in Eqs. (|3.7| a)-( |3-7| e) containing m s is infrared safe and therefore does not 
diverge in the limit m s — > 0. Also, one should take note of the fact that due to the 
equality of the heavy quark and heavy scalar diffusion coefficients, the more general 



, Q/ r = 32g*d A M 2 k 2 (l + cos 



+ 32g 4 d A M 2 k 2 (l - cos 



1 



(q 2 + m 2 ) 2 
1 

(q 2 + t7ij) 2 ' 



48g*d A M 2 k' 



(q 2 + m^) 2 



+Qg \ 2 = 8g 4 d A M 2 k 2 {l + cos 2 



+Qg \< = 8g 4 d A M 2 k 2 sm 2 6 



1 



(q 2 + m%f 



(q 2 + m 2 ) 2 ' 

| 2 = 8g*d A M 2 k 2 sm 2 6 1 2 ^ 

(q z + miV 



(3.7a) 
(3.7b) 
(3.7c) 
(3.7d) 
(3.7e) 
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heavy flavor diffusion coefficient D, given by the average of the two, coincides with 
Eq. ( |3.10|) as well. To better facilitate a comparison with the strong coupling limit of 
the theory — in which only the heavy flavor diffusion coefficient is a priori well-defined 
- we will in the following sections refer only to the latter quantity also in the weak 
coupling context. 

4. Discussion 

Having obtained an expression for the heavy quark (flavor) diffusion coefficient in 
weakly coupled Af = 4 SYM theory, it is interesting to analyze it and to compare 
it on one hand to the strong coupling result of Herzog et al. and others |Tl|, [12] and 
on the other hand to the corresponding weak coupling calculation in QCD by Moore 
and Teaney ||. An immediate observation from both the expanded result of Eq. ( |3.10| ) 
and the integral of Eq. ( |3.2| ) is that when written in terms of the 't Hooft coupling 
A = g 2 N c , the only explicit dependence on N c in the results comes from an overall 
factor of (1 — l/Af?) -1 . Keeping in mind that this multiplicative factor can be reintro- 
duced to the results at any later time, we shall in what follows, unless explicitly stated 
otherwise, consider the large N c limit and set (1 — = 1. 

In order to inspect the domain of validity of the small m D /T and m s /T expan- 
sions, we plot in Fig. our result for 1 / (DT) obtained both with and without the 
expansion, with the curve for the latter originating from a numerical evaluation of the 
integral in Eq. ( |3.2| ), in which the weak coupling expressions for the screening masses 
given in Eqs. ( |3.8| ) and ( |3.9| ) are used. We observe that the two curves begin to differ 
significantly at A ~ 3/4, and that the expanded result for D in fact starts to diverge 
when A > 2. This unphysical behavior signals the breakdown of the leading order weak 
coupling expansion for D and consequently implies that higher order corrections must 
be taken into account to gain even qualitative information about the intermediate cou- 
pling regime. Not performing the small screening mass expansion amounts to doing a 
partial resummation of our results, where some higher order contributions are included, 
but others, such as those coming from additional processes like bremsstrahlung or from 
corrections to the scattering amplitudes or screening masses, are neglected. While it is 
not a priori obvious that this is sufficient to gain quantitative information about the 
behavior of D at larger coupling, it is clear from Fig. [2] that including these corrections 
improves the qualitative behavior of the diffusion coefficient, as the unphysical diver- 
gence of D at A ~ 3 is removed. In what follows, we will therefore not use the small 
screening mass expansion, but instead evaluate the integrals of Eq. (|3.2j ) numerically. 

In Fig. Q.a, we investigate the behavior of our result at stronger coupling by plot- 
ting 1/(DT) as a function of A, with the screening masses still given by their leading 
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1/(DT) 




Figure 2: Plots of our two weak coupling results for 1/(DT): the expanded version taken 
from Eq. ( |3.10| ) (lower curve) and the one obtained via a numerical evaluation of the integral 
in Eq. ( |3.2| ) (upper curve). The large N c limit has been taken here, and the weak coupling 
expressions for the screening masses have been used. 



order weak coupling expressions. As A increases, an increasingly important source of 
ambiguity in this plot comes from neglecting the NLO corrections to the screening 
masses, which are expected to be sizable already at A ~ l. 3 Therefore, one must ex- 
ercise caution in interpreting these results and should preferably only use them when 
A <C 10. In Fig. |3|.b, we have circumvented this problem by plotting DX 2 as a function 
of m D /T, with rn^/m 2 ^ fixed, but have now replaced it with an ambiguity related to the 
choice of these ratios at any A > 1. From this figure it is, however, evident that the 
diffusion coefficient is relatively insensitive to deviations of the ratio m^/m 2 from the 
leading order weak coupling result of \ . 

For reasons of comparison, we have included already in Fig. |3|.a a plot of the strong 
coupling result for 1 / (DT) as obtained from Refs. [[IT], [l^] , according to which 

D = ^ (41) 

at large values of A. 4 As can be seen from Fig. §.a, it is easy to find a smooth, monotonic 
interpolating function that has the correct limiting behavior at small and large A, but at 

3 At these couplings, the leading order weak coupling results for the screening masses yield values 
> T . For consistency, one should therefore use resummed propagators in their evaluation, which would 
result in important (but at present unknown) 0(g 3 ) correction terms. For a discussion of this topic 
in the case of QCD, see Ref. p3| . 

4 The region of validity of this result is currently somewhat unclear, but the authors of Ref. [[LlJ 
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Figure 3: Left: A plot of the weak (lower curve) and strong (upper curve) coupling results 
for \/(DT) in the large iV c limit. The weak coupling curve has been obtained via numerical 
integration of Eq. ( |3.2| ), with the screening masses given by their leading order perturbative 
expressions, while the strong coupling curve is taken from Eq. (f4.1|). Right: The value of 
the (integrated) weak coupling result for DTX 2 as a function of m D /T for m^/rn^ = 0, 1/4, 
1/2 and 1 (from bottom to top) in the large N c limit. 

intermediate values of the coupling there is a wide region where neither result offers an 
accurate quantitative estimate for D. At A ~ 20, the weak coupling extrapolation of D 
is seen to be roughly six times larger than the corresponding strong coupling prediction, 
which is not surprising as we in any case are far beyond the region of validity of the 
weak coupling result here. However, we note that upon comparing the forms of the 
weak and strong coupling curves at intermediate couplings, it appears likely that the 
strong coupling expression yields an underestimate for the diffusion coefficient in this 
region. 

To put our discussion on a more quantitative footing, let us recall that using their 
weak coupling result, Moore and Teaney estimated the heavy quark diffusion coefficient 
for Nf = 3 QCD to be Dq « 1/T at a s = 0.5, where we have adopted the convention 
of denoting QCD quantities with tildes. In order to convert our result for the Af = 4 
diffusion coefficient at a = 0.5 into at least a rough estimate for this quantity, we refer 
to the behavior of the ratio Dq/ D at weak (and equal) coupling, as shown in Fig. |j. 
There, D is obtained by numerically integrating Eq. ( |B.3| ), with the screening masses 
given by their weak coupling expressions at N c = 3, while Dq is evaluated from the 



point out that at the experimentally interesting couplings of A ~ 20, it is already expected to obtain 
sizable corrections. Therefore, we urge the reader to use this expression only for A 3> 20. 
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Figure 4: Equal coupling plots of Dq/D for N c 
to bottom) as functions of a s . 



3 and N f = 0, 1, 2, 3, 4 and 8 (from top 



equivalent QCD expressions ||. We see from this plot that at least at weak coupling, 
the ratio is a slowly decreasing function of a s , and assuming this trend to carry on to 
stronger couplings, we estimate an upper bound Dq/D < 3 at a s = a = 0.5. Keeping in 
mind that the strong coupling result of Eq. (|4.1| ) is likely to be an overestimate at these 
couplings (corresponding to A ~ 19), we on the other hand obtain D\ a=05 > 1/(7T), 
which translates into the rough estimate Dq ~ 3/(7T) ~ 1/(2T) at a s = 0.5. 

Finally, from the point of view of the above comparisons between M = 4 SYM 
theory and QCD, it is also of some interest to investigate the ratio of the diffusion 
coefficients of the two theories at very weak coupling to get some insight into the order 
of magnitude of this quantity. Setting the couplings of the two theories again equal, 
we have from Ref. |8fl 



~ 72tt f 2T 
N f 



+ 



2N r 



, 4T 1 
log — + - 

m n 2 



1e + 



1e 



C(2) 
C(2) 

cm 

C(2) 



(4.2) 



from which we obtain at asymptotically weak coupling (where the logs in Eqs. Q3.10 ) 
and ([4.2|) dominate over the constant terms) 



We observe that for all reasonable values of Nj, the QCD diffusion coefficient is con- 
siderably larger than that of TV = 4 SYM, which is mostly a reflection of the fact that 
there are more light degrees of freedom in weakly coupled M = 4 SYM theory for the 
heavy quark to scatter off of. For example, in the Nf = case of pure Yang-Mills 
theory (in which all light particles are in the adjoint representation), one has 16 light 
bosonic degrees of freedom, while M = 4 SYM contains 64 bosonic and 64 fermionic 
degrees of freedom. A straightforward analysis shows that each bosonic degree of free- 
dom contributes equally to the leading log in k q , while each fermionic degree of freedom 
contributes half as much, so that at asymptotically weak coupling the diffusion coeffi- 
cient of pure Yang-Mills theory should be 6 ^. 32 = 6 times bigger than that of Af = 4 
SYM, just as we observed above. Similar conclusions have been drawn also in Ref. ||24|| , 
where the authors compared the weak coupling results for the sheer viscosity in M = 4 
SYM and QCD. 



5. Conclusions and future directions 

In the paper at hand, we have investigated the diffusion of a heavy, non-relativistic 
thermal particle — either a quark or a scalar belonging to a fundamental M = 2 hyper- 
multiplet — immersed in M = 4 Super Yang-Mills plasma. We have derived a result 
for the heavy flavor diffusion coefficient that is valid to leading order in g and T/M, 



and compared it to the corresponding strong coupling results of Refs. |L1|, [L2| as well 
as to the weak coupling calculations of Ref. || in QCD. Our findings show that a naive 
extrapolation of the weak coupling result to intermediate couplings yields a relatively 
large disagreement with the strong coupling predictions, while in the weak coupling 
limit the heavy flavor diffusion coefficient in the SYM theory is considerably smaller 
than the corresponding QCD quantity. Based on our analysis, we have estimated the 
heavy quark diffusion coefficient in QCD to be roughly Dq ~ 1/(2T) at a s = 0.5. 

As is evident from the small number of weak coupling results available in M = 4 
SYM theory, especially in comparison with the strong coupling limit or with pertur- 
bative QCD, there is a lot of further work to be done that can provide the QCD 
community useful insights from the abundance of existing AdS / CFT calculations. The 
obvious next goal related to the present work — and one that that should be straight- 
forward to achieve — is to generalize the results of this paper to the case of the diffusion 
of a relativistic quark with ^yv > 1. This is work in progress. 
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A. The Lagrangian 

In this first Appendix, we aim to present a somewhat detailed derivation of the La- 
grangian of our theory, M = 4 Super Yang-Mills with a massive M = 2 hypermultiplet, 
following to a large extent the treatment of Ref. ||18|| . The field content of the M = 4 
theory consists of one gauge multiplet with components (A^, A, D) and three chiral 
multiplets %■> x' an d x" with components (0, ip, J 7 ), (<//, ip', JF') and (</>", ip", J 7 "), respec- 
tively, while the M = 2 sector contains two fundamental massive chiral multiplets Q' 
and Q" with components and (<&",*", F"). Here, A^ is an SU(N C ) gauge 

field, A, ip, ip', ip", and W are Majorana fermions, <f), <f)', <ft", and $" are com- 
plex scalars, and D, JF, T\ T'\ F' and F" so-called auxiliary fields. All fields in the 
M = 4 sector transform in the adjoint representation of the gauge group and are hence 
traceless, hermitian N c x N c matrices, while the Af = 2 fields are fundamental under 
SU(iV c ) and can therefore be viewed as iV c -component vectors in color space. 

If we fix the heavy particle masses to M, the superpotential of the theory reads 



f(X, X', X", Q', Q") = -iV2 Q" T X Q' + 2iV2 tr ( X \x', x"}) + MQ' T Q". (A.l) 



Integrating out the auxiliary fields and going through some straightforward algebra, we 
obtain the Lagrangian |T8| 5 



5 In doing this, we have identified and corrected several misprints in the original reference. These 
are: an extra second term on the third-to last row of Eq. (27.4.1), a missing e matrix between A and 
ip on the second row of Eq. (27.9.3), reversed indices m and n in the second term on the tenth row of 
Eq. (27.9.33) and several misprints in the //-dependent terms of Eq. (27.9.33). The last two rows of 
the latter equation should be replaced by 



£ — Cq + C\ + C2, 



(A.3) 



where 



V2 Im (t' A fJ,) mn <i>A(t>T<i>m ~ V2 Im ([J.tf A )mn<f>A<f>m<f>n- 



m 



(A.2) 
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Co = -tr {-F^F^ + 2D^Dy + ^ + X0X + 2D^D^' + 2P M 0" t PV' 

- M 2 &' ] <&" - -ty'lpy' - - 2MRe vj/'P + ^", (A.4) 
Ci/g = Imj -4v / 2trAP+[0 t ,'0] + 4V2tr^'P + [0, V"] + V^M"^ [0', V] 

- 4v / 2tr^'P+[0", V] + 4V2tr V'F+[0' t , A] + 4^2tr ^"F + [0" t , A] 

- 2^ ^"P+^$ / - 2^ $ //T ^P+^' + 2^ $' t AP + ^' - 2^ ^"P + A$"* 

- 2 v / 2M$" t T $" - 2^2^$'^$' - 2V2 ^"P+0^'} , (A.5) 
£ 2 /<? 2 = -2tr I [<f>, 0'] 1 2 - 2tr | [0, 0*] | 2 - 2tr | [0, 0"] | 2 - 2tr | [0, <j>'«] | 2 

- i |trt a {2[0',0't]-2[0"t,0"]+W t -^ / ^ ,/T }r-tr[0 t ) 0] 2 

- 2 |trt a {2[0', 0"] + $'$ ,/T }| 2 -$ ,t {0,0 t }$' -$" T {0,0 t }$"*. (A.6) 

Here P± = | (1 ± 75), t a are the generators of SU(N C ) and a sum over a is implied. 

The form of the above functions can be greatly simplified upon making the redefi- 
nitions (adopted partially from Ref. ||19| ) 

'01 = V> ^2 = A, ^3 = 04 = 

= P+^' + P_^", (A.7) 

01 = =^=(X 1 + iY 1 ), 

02 = 0' =-^(X 2 +^ 2 ), 

03 = 0" = -^(A 3 + ^3), (A.8) 

$! = $ 2 = $"*, (A.9) 

where AT P and Y p are hermitian scalar fields and to a Dirac spinor. It is a straightforward 
exercise to show that in terms of these variables Lq reads 

Co = -tr {5 V" + AIP A + (DX P ) 2 + (DY p ) 2 } 

- $|,(-P 2 + M 2 )$ n - uj(Ip + iW>, (A.10) 
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where a summation over repeated indices is implied. Using the Majorana condition of 
Eq. (|1.2j) , the general form of the first six terms in Eq. ( |A.5| ) can on the other hand be 
written as 

4V2Imtr(^P + [0 fc ,%]) = -i2V2tr (^P+[0 P , ipj) - ^P-[Vi, 4 P }) 

= -2y/2tx (^[Re^ipj] - ipa 5 [lm(/) p , ty] J , (A.ll) 
which implies that one may simplify their sum into 



tr \ iipiO^JXp, ijjj] - ipa^ffi: [Y p , ijj^ 



(A.12) 



Here, a p and j3 p are coefficient matrices that may be taken as antisymmetric as the ipi 
anticommute and whose components can easily be verified to be given by Eq. ( ^.5| ). 
The remaining terms in Eq. (|A.5|) are simple to translate into the new variables, leading 
to the final result 



Ct/g = tr | - iipiafj [X p , if/j] + ipasPfj [Y p , ipj] } - u (Yi - 175X1) u 
+ 2v / 2Im| - cuP + ^i$i - $|-0iP + w + $J-0 2 P+^ - ^P + ^ 2 $2 
- 2M$t Fi$ n . 



(A.13) 



Finally attacking £2, the terms in Eq. ( |A.6|) that are independent of $ n read 



- 2tr 

- tr 



>1,<P2\ 



- 2tr [0! 



an 2 



2tr 



>1,¥>3 



2tr 



tr[0Ui] 2 



4tr I [0 2 ,. 



Using the Jacobi identity for the cross terms on the second line, we have 

-2 tr [02, 2 ] [03, 4] = 2 tr I [0 2 , 4] f + 2 tr | [0 2 , 3 ] | 2 , 
so Eq. ( |A.14|) becomes 



(A. 14) 



(A.15) 



- 2tr [0 1; 2 ] -2tr [0i, 



an 2 



2tr |[0 2 ,0 3 ]| 2 -2tr 

1 2 

-2 tr (HX^Xb]) , 



^2,4} 



2tr 



tr 



>i,</>3 



' - 2tr 
tr \q> 2 . o 2 



>n4 3 ]\ 2 



tr 



4,03] 



(A.16) 



where we have defined xa = (X, Yi, X2, Y2, X3, Y3). As before, the remaining terms in 
Eq. ( |A.6| ) are easy to translate into the new variables, leading to the result 



-17- 



A// = -^tr (i[ X a,Xb\? + (-iT^n ([02,4] + lh,4]) ®n 

- 4Re (^[02,03]^) " \ |(-l) n $t* a $n| 2 " 2|$^l| 2 
-$U01,0 f l}$n, (A.17) 

where repeated indices are again summed over. 
B. Matrix elements and integrals 

In this Appendix, we will briefly review our evaluation of the necessary scattering 
amplitudes squared in the non-relativistic limit, as well as explain, how one can perform 
the integrals in Eq. ( |3.2|) analytically in the weak coupling limit. Our treatment is to 
a large extent parallel to that of Ref . f| . 

B.l Matrix elements 

As discussed in Section || the scattering amplitudes squared for the heavy fermions and 
scalars become identical in the non-relativisic limit, which we exploit by only computing 
the simpler scalar amplitudes in the Coulomb gauge. We denote the color, flavor and 
momentum of the initial and final light particles by a, m, k and b, n, k', respectively, 
and the color and momentum of the heavy particles by i, p and j, p'. The angle between 
k and k' is written as 9, the structure constants of the gauge group as f a bc and the 
propagators for the scalars and the gauge field as G p cd and D^. Because the plasma 
has no preferred color orientation, we adopt the convention of averaging over the color 
configurations of the initial heavy particle, while the colors of the light particles as well 
as the final heavy scalar are summed over. 

As a concrete example, consider the amplitude for the process Sf — > Sf. The 
total scattering amplitude for this process is the sum of the first and fourth diagrams 
of Fig. [IJb and is given by 

Msf^sj = (-^mn/«fc«(k)7^(k / ))l>^(g)(^(P + P , ),(td)i J -) 

+ (^L/a^(k')75^(k)) GUQ) ( - 2igM(t d ) ij ) . (B.l) 

Upon squaring this expression and summing over m and n, it becomes evident that the 
cross term will be proportional to tr/5 1 , which vanishes due to the antisymmetricity 
of the matrix. Furthermore, in the non-relativistic limit we have (P + P') v ~ 2M5 u o, 
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so after summing over the colors, flavors and spins of the light fermions as well as the 
colors of the final heavy scalar, we obtain the result of Eq. (|3"7F|a), 



\M Sf ^sf\ = 32< ? 4 d A MV(l • cos„ !!i2y2 

1 



(q 2 + m 2 

+ 32g i d A M 2 k 2 (l - cos 9) . = - . (B.2) 



(q 2 + m 2 ) 



2 

sJ 



The results quoted in Eqs. (|3.7|b )- (|3.7|e ) are obtained in a highly analogous fashion. 
B.2 Integrals 

The integrals appearing in Eq. (|3.2|) are of the same type as those encountered in the 



QCD case, and our treatment of them therefore follows that of Ref. M quite closely. 
We begin by eliminating the three-dimensional delta function through integration over 
k', then change variables from p' to q = p' — p, and finally perform the angular part 
of the q integral to get rid of the energy delta function. This yields 



2k ( e f3k 



^ = ^M 2 L dk J d ™\^^^} M «^'\ 



e pk 



with f3 = 1/T. The sums over the amplitudes squared can be performed using the 
results of Eqs. (|377|a)-(|3~77|e) , and writing the results out explicitly we get 



SJ' 

, a 2 



+ 32 S * dA MV { ^) WT - l? , (B.4) 



q 2 q 4 



\M Hh ^ H v\ 2 = 8g 4 d A M 2 k 2 U - 

6,6' ^ 

J2 „A 



k 2 4A; 4 J (q 2 + m 2 ) 2 



~ ^d A M 2 k 2 i^-4a) ^ , , (B.5) 



k 2 4A; 4 J (q 2 + m 2 ) 2 

where the relation cos 8 = 1 — q 2 /2k 2 has been applied. 

The integral over q in Eq. ( B.3|) can be performed analytically, resulting in a some- 
what lengthy one-dimensional integral representation for the diffusion coefficient as a 
function of the screening masses, which we plotted numerically in Section 4. In the 
true weak coupling limit, where the screening masses satisfy m/T ~ j < 1, we may, 
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however, simplify the calculation considerably by noting that all terms in Eqs. (|B.4|) 
and ( p.5|) proportional to positive powers of q are infrared insensitive and thus inde- 
pendent of the masses to leading order. This enables us to set the masses to zero in 
these terms and gives 

E \MH f -+Hf> I 2 = eVd^fc 2 1 , (B.6) 
/,/' yq d) 

£ |AW,| 2 = S 9 'd A MV (^L _ _L + , (B.7) 

which simplifies the result of the q integration in Eq. (|B . 3|) dramatically. Finally eval- 
uating the k integrals with standard methods, our weak coupling result for 3k h reads 

g A d A f°° JIl2 f e? k ( , 4£; 2 \ e? k ( 3 , 4& 2 \ 1 

3k h = / dfcfc 2 4 ? -1 + log — + ? — + log — ^ 

2tt 3 y \( e /3fe + l) 2 V (e^-l) 2 V 4 m 2 ; J 

3 4 d A T 3 f , 2T 13 1 C'(2) 1 ,„ , 

J log— +— - 7£ + -log2 + ^-f L (B.8) 



2vr \ m D 12 m 3 ta C(2) J ' 
which in turn leads to the expression of Eq. ( |3.10| ) for Dq (or D$ 
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